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ABSTRACT 

HERSCHEL-HIFI observations of water from the intermediate mass protostar NGC7129 FIRS 2 provide a powerful diagnostic of the physical 
conditions in this star formation environment. Six spectral settings, covering four H2 16 and two H2 18 lines, were observed and all but one H2 18 
line were detected. The four H 2 16 lines discussed here share a similar morphology: a narrower, ss 6 km s _I , component centered slightly redward 
of the systemic velocity of NGC7 129 FIRS 2 and a much broader, sa 25 km s component centered blueward and likely associated with powerful 
outflows. The narrower components are consistent with emission from water arising in the envelope around the intermediate mass protostar, and 
the abundance of H2O is constrained to «10~ 7 for the outer envelope. Additionally, the presence of a narrow self-absorption component for the 
lowest energy lines is likely due to self-absorption from colder water in the outer envelope. The broader component, where the H 2 0/CO relative 
abundance is found to be » 0.2, appears to be tracing the same energetic region that produces strong CO emission at high J. 
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1. Introduction 

Observations of star formation in nearby molecular clouds have 
provided a reasonably clear picture of this process for low-mass 
stars. Nevertheless fundamental questions remain. One of these 
is the variation in the formation process between low through 
high mass stars. The intermediate mass (IM) sub-program within 
the Water In Star-forming regions with HERSCHEL (WISH) 
key program (van Dishoeck et al 2010) aims to use observations 
of water, and complementary molecules, to probe the regions 
around several intermediate mass protostars - focusing on the 
structure of the envelope and the energetics of the outflow. These 
sources are usually defined as young stellar objects (YSOs) with 
bolometric luminosities between 75 L and 2 x 10 3 L . A key 
goal for the longer term of this investigation is to place IM star 
formation in context with its low- and high-mass brethren. 

An excellent example of an extremely young IM protostar is 
NGC7129 FIRS 2 (Eiroa et al 1998; Fuente et al EOTJll 120051 
located at a distance of 1260 + 50 pc from the Sun (Shevchenko 
& Yakubov, |1989l l. With a luminosity of 500 L Q and estimated 
stellar mass of 5 M Q , FIRS 2 lies near the middle of the IM 
luminosity range. The protostar has produced a hot core (Fuente 
et al[2005 ) but has no evidence for a large, well-developed disk - 
implying it is a young source. A powerful quadrupolar outflow, 



* HERSCHEL is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im- 
portant participation from NASA. 



likely due to the superposition of two bipolar jets, is found very 
close to this source (Fuente et al 2 0011 >. 



NGC7129 FIRS 2 has been the subject of significant inves- 
tigation over the last year. A robust model for the enshroud- 
ing, and pre-natal, envelope has been proposed by Crimier et 
al (2010 ) utilizing all available far infrared (Spitzer) and submil- 
limeter (JCMT, IRAM) brightness measurements. Additionally, 
Fich et al (2010) used HERSCHEL PACS observations of highly 
energetic CO (up to CO J = 33-32 which arises from a level at 
an energy equivalent 3093 K) to analyse the energetics associ- 
ated with NGC7129 FIRS 2. They found that these CO lines are 
much brighter than expected from the envelope alone, revealing 
the need for additional heating beyond simple reprocessing of 
the protostellar radiation. A warm slab model with temperatures 
greater than 1000K, considered a proxy for shock heating in or 
along the surfaces of the outflow lobes, provided a much bet- 
ter fit to the CO observations. A lack of spectral resolution in 
the PACS data hindered further investigation into the connection 
between the CO lines and the outflow. 



In this paper we present a first-look at the water spectra 
obtained toward NGC7 129 FIRS 2 with the HIFI insttument 
(de Graauw et all2TJT0l ) on the HERSCHEL Space Observatory 
(Pilbratt et al 2010). The wealth of spectrally resolved water 
lines observed provide powerful diagnostic measures for the en- 
velope and the outflow. 
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Fig.l. H 2 spectra detected toward NGC7129 FIRS 2 with 
HIFI. Note the similarities in morphology between the lines, es- 
pecially the appearance of at least two emission components, 
narrow and broad (see Figure [3]). The lowest energy transi- 
tion lines also show evidence for somewhat blue-shifted self- 
absorption. 



NGC7129 FIRS 2, with the most likely source of absorption be- 
ing the protostellar envelope. Two CO isotopologues were also 
detected (see Figure |2). While both these lines peak at the sys- 
temic velocity and are much narrower than the H2 16 lines, the 
13 CO 10-9 line is broader than the C ls O 5-4 line (see Table|2]i. 
Of the water isotopologues, only the H2 ls O lio - loi line was 
sufficiently strong to be unambiguously detected (Figure|2]l. 

Table 2. Water and CO line emission toward NGC7129 FIRS 2 



2. Observations 

HERSCHEL HIFI was used to obtain data for this project dur- 
ing the HERSCHEL science demonstration phase and the HIFI 
priority science program, from 3-20 March, 2010. The obser- 
vations were taken using the fast dual-beamswitch mode with a 
standard nod of 3' using receiver bands 1, 2, and 4. In total six 
spectral settings were selected, covering four H2 le O transitions 
and two H2 ls O transitions, with integration times between 15 
minutes and 1 hour. Two isotopologue lines of CO were also 
located within these settings. CO 10-9 was also observed on 
June 15, 2010, with a 10 minute integration in band 5 using 
the fast dual-beamswitch mode. The complete set of observed 
lines are shown in Table Q] All the observations were taken 
toward NGC7129 FIRS2, at 21 /, 43'"1.7' s 66°03'23" (J2000). 
Diffraction-limited beamsizes range from 20" - 40" (equivalent 
to 25,000 - 50,000 AU at the distance of NGC7129). 

Table 1. Observed H 2 le O, H^O, and CO transitions. 



Transition 


V 


EJk B 


Transition 


V 


EJk B 




(GHz) 


(K) 




(GHz) 


(K) 


p-H 2 "'0 ln-Ooo 


1113.34 


53.4 


o-H 2 ls O l 10 -loi 


547.68 


60.5 


p-H 2 16 O 2 02 -l„ 


987.93 


100.8 


p-H 2 18 O 2 02 -l n 


994.68 


100.7 


p-H 2 16 O2 n -2 02 


752.03 


136.9 


C ls O 5-4 


548.83 


79.02 


o-H 2 16 O3 12 -3 3 


1097.37 


249.4 


13 CO 10-9 


1101.35 


290.79 








CO 10-9" 10-9 


1151.99 


304.2 



<a) Only a preliminary reduction of this observation has been performed, it is discussed only 
in Section 3.2 

The data were pipelined using HIPE 3.(fl (Ott 2010) and the 
data analysis was completed using both HIPE 3.0 and CLASS. 
The results of the two analyses were found to be the same to 
better than twenty percent, where the largest uncertainties were 
found during the fitting of multiple components. Before the re- 
duction, the chop positions were differenced to check for ex- 
tended emission and none was detected. During the reduction, 
the H and V polarizations were compared to assess quality and 
then averaged. In most cases the two polarizations showed no 
significant variations. The H2 lfi O 3 12— 3o3 line, however, showed 
both a lower peak intensity and line width for the V polarization, 
resulting in an 8% difference in integrated intensity. For both po- 
larizations, the overall morphology of the spectrum is similar, 
with the variation likely due to the extended nature of the source 
and slight offsets in the H and V beam positions. Finally, the 
WBS and HRS spectra were compared for consistency. Regions 
containing spurs were excluded from the analysis and a first or- 
der polynomial baseline was removed from each spectrum. The 
WBS data were ultimately smoothed to the goal resolution of 
1 . 1 MHz and the antenna temperature T* A was converted to T m \, 
by applying a main beam efficiency of 0.74. 

The four H2 ls O transitions were observed and are shown in 
Figure Q] The spectra show evidence for at least two compo- 
nents: a narrower line peaked slightly redward of the systemic 
velocity of the source (-9.8 km s _1 ) and a broad line blue-shifted 
from the systemic velocity. The ground-state transition shows 
significant self-absorption, with a minimum at -10.4kms _1 , 
hints of which can be seen in the other observed transitions, and 
is likely due to colder water lying along the line of sight toward 



1 HIPE is a joint development by the HERSCHEL Science Ground 
Segment Consortium, consisting of ESA, the NASA HERSCHEL 
Science Center, and the HIFI, PACS and SPIRE consortia. 



Transition Component rms° T^ b [ T MB do '' Vlsr 1 ' FWHM 7 







(mK) 


(K) 


(Kkmr 1 ) 


(kmr 1 ) 


(km s" 


H 2 16 Iu-Oqo 


no Gaussian fit 


16 


0.59 


8.4 


-8.2 




2o 2 -ln 


naiTow 


20 


0.36 


2.3 


-7.2 


6.2 




broad 


20 


0.23 


6.1 


-12.4 


24.6 


2n-2o 2 


narrow 


16 


0.21 


1.5 


-7.5 


6.5 




broad 


16 


0.11 


3.1 


-13.7 


26.8 


3i 2 -3o3 


narrow 


16 


0.14 


1.7 


-7.3 


7.8 




broad 


16 


0.11 


2.8 


-14.5 


24.4 


H 2 ls O lio-loi 




5 


0.015 


0.071 


-8.5 


5.3 


H 2 IS 2„ 2 -1„ 


not seen 


5 




<0.010' 






C is O 5-4 




5 


0.16 


0.49 


-9.8 


2.1 


13 CO 10-9 




20 


0.13 


0.44 


-9.8 


5.5 



(<Y| The rms measure is taken over a 1 km s -1 range. 

^ Uncertainties in the peak temperature and integrated line intensity are dominated by 
calibration issues. For H? ]6 Iu-Oqo, 2q2— In. 3i2-3o3, and 13 CO 10-9 the uncertainty is 
32%. For all other lines the uncertainty is 16%. 
' c) Velocity accuracy better than -0.5 km s -1 . 

A 20% uncertainty in the line widths is due primarily to the component separation. 
' e) The 3o" detection limit, for data smoothed to the expected line width (5 km s -1 ), is 0.015 K 



3. Discussion 

The H 2 16 lines detected toward NGC7129 FIRS 2 provide an 
opportunity to examine the physical and chemical environment 
of this young IM protostar. Attempts to fit the spectra with two 
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Fig. 2. Isotopologue CO and water spectra observed toward 
NGC7129 FIRS 2 with HIFI. Note the width of the isotopo- 
logues relative to the main line. 
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Fig. 3. Spectrum of the H2 16 2o2 - In transition (blue) su- 
perposed with the two Gaussian fit (black) from Table [2] The 
fit includes a narrow, 6.2 kms" 1 FWHM, component centered 
at -7.2 kms" 1 (grey) and a much broader, 24.6 kms" 1 FWHM 
component, shifted slightly blue-ward at -12.4km s" 1 (purple). 



Gaussian components proved reasonably successful for the ex- 
cited lines (Table [2] & Figure |3j. For these three lines, the nar- 
rower component was well described by a * 6 kms" 1 , Gaussian 
centered redward of the systemic velocity of the source near 
-7 kms" 1 coupled with a broader, 25 kms , Gaussian cen- 
tered near -13.5 kms 1 (Figures [TJ& EJ- The likely additional 
absorption component in each of these spectra could account 
for the small variations in the line centroids and widths of the 
fits. Similar profiles are seen for low- (Kristensen et al. 2010) 
and high-mass (Chavarna et al. 2010) protostars, where they 
are named medium and broad components. Even narrower lines, 
« 1-2 kms -1 FWHM, are seen in self-absorption of the lowest- 
energy water lines as well as in ground-based emission lines 
such as CS 3-2. The latter profile also reveals a 6 km s" 1 compo- 
nent (Fuente et al. 2005). 



Beyond the morphological nature of these lines, it is possible 
to compare the water emission against model predictions. Fich et 
al ( 1201 Oi l note that NGC7129 FIRS 2 has two main components 
which may account for the observed emission. The prenatal en- 
velope, which is warmed by the IM protostar embedded within 
it, may produce water emission. Also, the source powers an en- 
ergetic outflow which could produce water emission either in the 
outflow itself or along the shock heated walls, where the outflow 
interacts with the denser envelope. In the following subsections, 
we consider the water emission from each of these components 
and compare the expectations against the HIFI observations. 

3.1. Water from a spherical protostellar envelope 

The spherical envelope model for NGC7129 FIRS 2 computed 
by Crimier et al (2010) fits the spectral energy distribution from 
the far infrared through the submillimeter very well and thus pro- 
vides a reasonable set of physical conditions for examining the 
observed water emission. This model has an envelope mass of 50 
M Q , an optical depth at 100 fim of 2.3, an inner radius of 100 AU 
and an outer radius of 18,600 AU. The temperature at the inner 
envelope radius is 289 K, falling to 100 K at a radius of 373 AU 
where the H2 density is 4.4 x 10 7 cm -3 . The density varies as a 
power-law with index -1.4. 

Model water lines were calculated using the RATRAN code 
(Hogerheijde & van der Tak 2000) using collisional cross sec- 
tions for H20-He from Green (1993), scaled by 1.348 to make 
a first order approximation for collisions with H2. An ortho-to- 
para ratio of 3 was assumed for H2. The density and temper- 
ature profile of the Crimier et al ( 120101 1 model were followed 
and the radial velocity profile was assumed to be the free-fall 
velocity appropriate to a central object of mass 1 . 1 M Q (ie. ve- 
locity of -4.2 kms" 1 at the inner radius and -0.31 kms" 1 at the 
outer edge). Models were calculated for two abundance zones 
in the envelope: a warm, > 100K inner envelope with a higher 
abundance and a cooler outer envelope with a lower abundance, 
perhaps characteristic of a "freeze-out" zone. 

As a starting point, we modelled the ortho-H2 ls O lio - lot 
observation. This line has a similar profile to the narrower com- 
ponents of the H2O lines and if it can be reproduced by the spher- 
ical envelope model we may conclude that the narrower compo- 
nent of the H2O lines also are likely to arise in the envelope. A 
large number of models were run over a range of inner (10" 9 - 
10" 4 ) and outer (10" 7 -3xl0" n ) envelope abundances. The non- 
thermal velocity width parameter b was also varied and values 
> 3 kms" 1 were strongly ruled out, while b = lkms" 1 pro- 
duced too narrow a line. The best fit was found for b — 2 km s" 1 . 

The abundance of ortho-H2 18 in the outer envelope is well- 
constrained to be 3 + 1 x 10" 10 , while inner envelope abun- 
dances between 3 x 10" 7 and 1 x 10" 5 produce results consis- 
tent with the integrated flux of the ortho-H2 18 lio - lot line 
(Table IaTT ). These are surprisingly high values; indeed, the top 
end of this range is impossible when the cosmic ls O abundance 
is considered. The para-H2 1 s O 2 2 - 1 1 1 was modeled in a similar 
way. The non-detection of this line constrains the outer envelope 
abundance of para-H2 18 to < 3 x 10" 11 . The model fits are in- 
sensitive to the inner envelope abundance used (Table lA~2l ) and 
this may be a consequence of beam dilution. 

Assuming that the ratio of 16 0/ 18 = 550 is also representa- 
tive of the ratios of both of the ortho and para H20/H2 I8 0, we 
examined models for the narrower components of the H2O lines 
with outer envelope abundances near 1 x 10" 7 and inner abun- 
dances up to 3 x 10" 4 (limited by the cosmic O abundance). The 
best fitting models indicate that the outer envelope abundance of 
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total H2O is of order a few 10" 7 , and are summarised in Table 
[3] The ortho to para ratio implied by the best fit outer envelope 
abundances is 10:1, this is highly unlikely and within the calibra- 
tion (<32%) and model uncertainties the ratio is consistent with 
3. The abundance of the inner envelope could not be constrained 
and the observed line profiles were poorly reproduced, particu- 
larly in the case of the para lines. Whilst the parameter space of 
the models are still not fully explored, this suggests limitations 
in the model used with the assumptions concerning the velocity 
field and spherical structure most suspect. 

Table 3. Overview of best fitting spherical envelope models 



Isotopologue 


Outer abundance a 


Comment 


o-H, ls O 


3 ± 1 X icr lu 


implies o-H 2 1.7 ± 0.5 X KT' " 


p-H 2 ls O 


<3xl(T n 


implies p-H 2 1.7 X 1(T 8 * 


o-H 2 


3 x 1(T 7 


Well constrained by both excited para lines 


p-H 2 


3 x 1(T 8 


> 10 -8 required for self-absorption 



(H) The relative abundance of species listed, and not total H2O 
(4) Assuming 16 0/ ls O = 550 



3.2. Water from a slab model 

Fich et al 120 101 concluded that the energetic PACS CO lines ob- 
served toward NGC7129 FIRS 2 could not be fit by a reprocess- 
ing spherical envelope model as there was an insufficient volume 
of warm gas to produce the strong emission observed. It was 
postulated that the CO emission might arise along shock heated 
walls in the outflow cavity, where the gas temperature can reach 
extreme values. NGC7129 FIRS 2 is known to contain a pow- 
erful outflow (Fuente et al 2001) and observations of low-mass 
stars have shown that such outflows can produce heated walls 
via UV photons (van Kempen et al 2009) and shocks (Gianninni 
et al 1999, van den Ancker et al 2000, Nisini e t al 20 02, Arce 
et al 2007, van Kempen et al 2010). Fich et al d20T0b used the 
RADEX code to show that the observed CO emission could be 
explained within a model of a slab geometry with a temperature 
of 1100 K, an H 2 density of 1.0 x 10 8 cm 4 , and a CO column 
density of 10 14 cirT 2 . This model constrained the lowest temper- 
ature at which the CO observations could be adequately fit. A 
more comprehensive set of CO models was produced for this 
paper, and a best fit is found for a temperature of 1200 K, an H2 
density of 10 7 cirT 3 and a column density of CO of 1.7 x 10 14 
cm -2 . Assuming a CO abundance of 10" 4 , these results suggest 
a slab with thickness dz ~ 10 10 cm. 

Using these same conditions (1200 K and 10 7 cm" 3 ), but 
substituting water for CO, may provide an additional constraint 
on the possible shock heating conditions, and the conditions un- 
der which the broad component of the water lines arise. The wa- 
ter lines observed with HIFI set the limit on the ortho-H20 col- 
umn at 2.3 + 0.2 x 10 13 cm" 2 while the para-H20 lines observed 
set a maximum limit in the column at 7.5 + 0.5 x 10 12 cm" 2 , 
giving an ortho-to-para ratio of 3. 

These models imply a water abundance that is approximately 
five times lower than CO and can be compared with that found 
from analysis of the H2O and CO line ratios in the line wing, 
where the emission is optically thin. The CO 10-9 line exhibits 
a similar profile to that seen in the H2O lines (Figure IA.U and 
be compared directly with the H2O 2(>2 - 1 1 1 line, as they are 
observed in almost the same beam (22" vs 19"). The abundance 
ratio was calculated with RADEX using a density of 10 6 cm" 3 . 



Above T- 150 K the results are insensitive to changes in temper- 
ature. The abundance ratio is found to be « 0.2 near the line cen- 
tre and approaches 0.3 at the highest velocities in the red wing 
(Table lA~3l ). This is consistent with the findings from the slab 
model and also with the behaviour seen in the low-mass sources 
of Rristensen et al. (2010), who report this indicates that « 10% 
of the available O is in H2O. 

That the same model can be used to fit both the PACS CO 
observations and the broad component of the HIFI H2O lines 
supports the hypothesis that these lines arise along shock heated 
walls in the outflow cavity. Shock models will allow to better 
constrain the conditions under which these broad lines arise. 

4. Conclusions 

HERSCHEL HIFI spectroscopy of water in the vicinity of the 
intermediate mass protostar NGC7129 FIRS 2 has revealed that 
the water emission arises from at least two sources. The ob- 
served emission lines can be decomposed into both a narrower, w 
6kms"', component and a much broader, « 25kms" 1 compo- 
nent. The integrated intensity of the narrower component can be 
fit by a simple free-falling envelope model with an outer enve- 
lope total H2O abundance of « 10" 7 , although the shape of the 
line profiles are not reproduced. The broader component appears 
to be related to the heated gas already observed in high-/ CO 
lines with PACS and associated with the known energetic out- 
flow. These initial modeling results suggest that a dedicated pa- 
rameter study for this source should prove extremely fruitful in 
constraining the physical and chemical conditions in NGC7129 
FIRS 2. 
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Appendix A: Online Material 

Table A.l. Summary of o-H 2 18 1 io - 1 oi model results. 



Envelope abundance 
Inner Outer 


ypeak 
MB 

(K) 


/ Tun dv 
(Kkms- 1 ) 


lxlCT 3 


3x10-'" 


0.017 


0.071 


3x10-" 


3xl(r 10 


0.018 


0.075 


1x10-" 


3x10-'° 


0.017 


0.069 


3xl0" 7 


3x10-'° 


0.016 


0.065 


lxicr 7 


3x10-'° 


0.015 


0.060 


3xl0" s 


3x10-'° 


0.014 


0.054 


lxl(T 8 


3xl(r'° 


0.014 


0.050 


lxiir 3 


2x10-'" 


0.015 


0.063 


3xl0" 6 


2x10-'° 


0.014 


0.060 


lxlO" 6 


2x10-'° 


0.014 


0.060 


3xl0- 7 


2x10-'° 


0.013 


0.052 


lxiir 3 


1x10-'" 


0.010 


0.043 


1x10-" 


lxKT 10 


0.009 


0.038 


3x10-' 


1x10-'° 


0.009 


0.036 


lxl(T 7 


1x10-'° 


0.008 


0.032 


3xl0- 8 


1x10-'° 


0.007 


0.027 



(0) The observed line has a 7* e jf = 0.015 K and an / T MB dv = 0.071 KknT 1 , the best fitting 
model is selected by comparison with the latter. 
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Fig. A.l. The H2O 2 ( )2 -In and CO 10-9 transitions; the main 
isotopologue of CO exhibits a profile similar to the main iso- 
topologues of H 2 0. 



Table A.2. Summary of p-H 2 O 2o2 - In model results, which 
are compared with J T M b dy < 0.010 K km -1 



Envelope abundance 


.j. peak 


/ T M b dv 


Inner 


Outer 


(K) 


(Kkms-') 


3x10-" 


1x10-'" 


0.0088 


0.033 


lxlO" 6 


1x10-'° 


0.0088 


0.033 


3xl0~ 7 


1x10-'° 


0.0085 


0.032 


3x10-" 


3x10-" 


0.0026 


0.0097 


lxlO"" 


3x10-" 


0.0027 


0.0099 


3xl0- 7 


3x10-" 


0.0023 


0.0085 


3x10-" 


1x10-" 


0.0008 


0.0032 


1x10-" 


1x10-" 


0.0008 


0.0031 


3xl0- 7 


1x10-" 


0.0008 


0.0030 



Table A.3. H 2 2 02 - 1 1 i/CO 10-9 abundance ratios in 5 km s" 1 
intervals, calculated from CO 10-9/ H 2 2 2 - In line ratios, 
using n = 10 6 cm" 3 and T > 150 K 



dv L SR 






(kms- 1 ) 


CO 10-9/ 


x(H 2 0)/ 




H2O 2o2 - 111 


x(CO) 


-10 --5 


3.25 


0.27 


-5-0 


3.86 


0.22 


0-5 


4.14 


0.21 


5-10 


1.86 


0.47 



